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1.   INTRODUCTION 

In  recent  years,  the  theory  of  networks  has  received  in- 
creased interest  in  the  field  of  operations  research,  in  con- 
nection with  the  solution  of  allocation  problems,  shortest  route 
problems,  and  other  flow  problems.   T.  C.  Koopmans  (9),  in  his 
work  on  transportation  problems,  was  the  first  to  interpret 
properties  of  optimal  and  nonoptimal  solutions  with  respect  to 
the  linear  graph  associated  with  a  network  of  routes.   This  work 
was  then  followed  by  the  pioneering  work  on  flows  in  networks 
by  Ford  and  Fulkerson  (6,  7) . 

The  purpose  of  this  report  is  to  present  the  application 
of  the  theory  of  networks  to  transportation,  transshipment,  and 
maximal  flow  problems.   In  Section  2  some  definitions  and 
theorems  concerning  the  theory  of  network  and  maximal  flow  prob- 
lems are  presented. 

In  Section  3  a  simple  transportation  problem  is  solved  by 
four  methods.   The  type  of  transportation  problem  considered 
consists  of  a  number  of  source  points  where  the  commodities  are 
available  in  limited  quantities,  and  a  number  of  demand  poinds 
where  there  is  a  requirement  for  the  commodity.   The  cost  of 
shipping  one  unit  of  commodity  from  a  source  point  to  a  demand 
point  is  known.   The  problem  is  to  determine  the  minimum  cos; 
route  of  transporting  the  commodity  and  to  satisfy  all  the  re- 
quirements.  The  first  method  used  to  solve  the  example  is  the 
simplex  algorithm,  developed  by  Dantzig  (2).   The  second  is  the 
one  commonly  called  "the  transportation  algorithm"  (2,  8,  10), 


which  is  a  modification  of  the  simplex  method.   The  third  method 
is  the  discrete  maximum  principle  which  has  been  shown  by  Pan, 
Wang,  and  others  (1,  $) ,    that  it  is  applicable  to  solving  trans- 
portation problems.   The  fourth  and  final  method  used  in  sec- 
tion 3  is  network  theory  as  applied  to  the  transportation  problem. 

In  section  I4.  a  transshipment  problem  is  solved  to  point 
out  the  ease  of  solving  problems  with  networks  as  compared  to 
the  simplex  method  and  the  transshipment  algorithm  (2).   The 
transshipment  problem  is  a  special  type  of  transportation  prob- 
lem where  several  possible  routes  for  shipping  exist  between  a 
source  point  and  a  demand  point.   The  main  difference  and  the 
added  difficulty  of  transshipment  problems  is  that  these  inter- 
mediate points,  where  neither  supply  nor  any  demand  exists,  have 
to  be  considered.   An  example  problem,  consisting  of  eight  nodes, 
two  source  points,  three  demand  points,  and  three  intermediate 
points,  is  set  up  as  a  linear  programming  problem,  and  is  solved 
by  the  transshipment  algorithm  and  the  method  of  networks.   All 
the  methods  used  in  this  report  are  applicable  only  to  the 
problems  with  linear  cost  function,  except  discrete  maximum 
principle . 

The  application  of  the  theory  of  maximal  flows  in  networks 
is  presented  in  section  $.      The  maximal  flow  problem  is  one  of 
determining  the  maximum  flow  from  a  source  to  a  sink  when  each 
arc  of  the  network  has  some  fixed  flow  capacity.   The  sink  is 
the  ultimate  destination  of  flow  originating  at  the  source. 
The  commodity  which  flows  in  networks  may  be  such  items  as 
fluids,  electricity,  funds,  and  automobiles,  where  the  quality 


is  measured  as  rate.   An  example  is  solved  using  the  method  of 
chains  of  positive  arc  capacities,  the  labeling  method,  and 
the  matrix  method.   These  methods  are  presented  in  Ford  and 
Pulkerson  (6,  7)  and  basically  depend  on  the  max- flow  rain-cut 
theorem  (3,  k)  ,    which  apply  to  steady-state  flows,  that  is, 
arc  flows  which  do  not  change  with  time.   It  is  noted  that 
dynamic  programming  and  discrete  maximum  principle  can  be 
applied  to  all  the  problems  considered  in  this  report  but  not 
much  work  has  been  done  at  this  time. 


2.   SOME  DEFINITIONS  AND  THEOREMS 

Network.   A  network  or  linear  graph  consists  of  a  number 
of  nodes,  junction  points,  or  points,  each  joined  to  some  or 
all  of  the  others  by  arcs,  links,  branches,  or  edges. 

The  numbers  with  circles  in  Pig.  1  are  the  nodes.   The 
arcs  or  branches  are  shown  by  straight  or  curved  line  segments, 
each  of  which  links  two  nodes,  e.g.,  (1,  2),  (1,  3),  (3,  2), 
(5,  h,),    etc. 

If  i  and  j  are  two  nodes,  then  (i  —  j)  is  used  to  denote 
a  directed  arc  and  represents  an  allowable  precedence  between 
i  and  j.   In  Fig.  1,  if  it  is  allowed  to  proceed,  starting  at 
1,  to  3  or  k;    °r  starting  at  2,  to  1  or  3;  or  starting  at  3,  to 
2  or  5;  or  starting  at  k,    to  $,    then  the  directed  arcs  are 
(1  —  3),  (1  -»  k),    (2  —  1)»  (2  -*  3)>  (3  -*"  2)'  (3  -  5),  and 
(1;  _>  5),  which  form  a  directed  or  oriented  network  as  shown  in 
Fig.  2.   Arcs  in  the  case  of  a  directed  network  are  also  known 
as  oriented  arcs  or  oriented  branches. 

Chain.   A  chain  is  a  sequence  of  arcs  (i,  ix) ,  (ii,  i2), 
.,  (ivi  j),  connecting  the  nodes,  i  and  j,  regardless  of 
the  ways  in  which  these  arcs  may  be  directed  (Fig.  3). 

Path.   When  every  arc  of  a  chain  has  the  same  direction, 
the  chain  so  formed  is  called  a  path  (Fig.  l±)  . 

Loop.   A  chain  of  arcs  connecting  node  i  to  itself  is 
called  a  loop  (Fig.  5) • 

Tree.   A  network  having  no  loops  in  which  every  node  is 
connected  to  every  other  node  through  a  chain  of  arcs  is 


Fig.  1.   Example  of  a  network  (undirected) 


Pig.  2.   Example  of  a  directed  network. 
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Pig.  3-   Example  of  a  chain. 
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Pig.  1|.   Example  of  a  path. 
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Fig.  5.   Examples  of  loops. 


Fig.  6.   Example  of  a  tree  (directed) 


called  a  tree  (Pig.  6) . 

Theorem  1.   A  network  having  n  nodes  is  a  tree  if  it  has 
(n  -  1)  arcs  and  no  loops. 

Capacity.   Capacity  of  the  arc  (i  ->  j )  is  the  maximal 
amount  of  some  commodity  that  can  arrive  at  node  J  from  node  i 
per  unit  time.   Hote  that  it  is  possible  to  have  infinite  arc 
capacities. 

Flow.   Plow  in  the  arc  (i  — >■   j)  is  the  amount  of  some  com- 
modity allowed  to  arrive  at  node  j  from  node  i  per  unit  time 
and  is  always  less  than  or  equal  to  the  capacity  of  that  arc. 
Directed  or  oriented  arcs  represent  the  direction  of  flow  in 
that  arc. 

Source,  Sink,  and  Intermediate  Nodes  or  Points.   A  node  i 
is  called  a  source  if  every  arc  which  has  this  node  as  an  end 
point  is  oriented  in  such  a  way  that  the  flow  in  the  arc  moves 
away  from  i  to  another  node.   A  node  j  is  called  a  sink  if 
every  arc  which  has  this  node  as  an  end  point  is  oriented  in 
such  a  way  that  the  flow  is  from  other  nodes  to  j.   And  all 
other  nodes  are  known  as  intermediate  nodes. 


3.   TRANSPORTATION  PROBLEM 

A  classical  transportation  problem  can  be  described  as 
follows.   A  given  quantity  of  a  commodity  is  available  at  each 
of  several  origins  or  sources.   There  is  a  requirement  for  the 
commodity  at  each  of  the  destinations  or  sinks.   The  cost  of 
shipping  one  unit  from  an  origin  to  a  destination  is  known. 
Also  it  is  possible  to  ship  from  any  origin  to  any  destination. 
The  costs  of  shipping  between  the  sources  and  the  destinations 
are  linear  and  proportional  to  the  amount  shipped.   The  problem 
is  to  determine  the  minimum  cost  routing  from  the  origins  to 
the  destinations  and  to  satisfy  all  the  requirements. 

The  transportation  problem  can  be  put  in  equation  form 
as  follows: 

Let  there  be  m  origins  (supply  points)  and  n  destinations 
(demand  points) . 

Let  iij  be  the  number  of  units  shipped  from  origin  i  to 
destination  j;  let  0^   be  the  number  of  units  of  the  commodity 
available  at  origin  i,  and  D_.  the  number  of  units  required  at 
destination  j.   Let  Cjj  be  the  cost  of  shipping  one  unit  from 
origin  i  to  destination  j.   The  standard  transportation  algo- 
rithm is  based  on  the  assumption  that  the  total  demand  equals 
the  total  supply.   Of  course,  it  is  always  possible  to  satisfy 
this  condition  by  inserting  dummy  demand  or  supply.   Thus  the 
problem  is  to  find  Xj  •  J:  0  which  minimizes  the  objective 
function, 


m    n 

1=1  j=i 

subject  to  the  supply  constraints, 


(1) 


33 

n. 


=  Qt,  i  =  1,    2,     .     .     .,    m        (2) 


3=1   1J 

and  the  demand  constraints, 

H    x.  .  =  D.,       j  =  1,  2,  .  .  .,  n      (3) 

1=1   J    J 

with  the  condition  that  the  total  demand  equals  the  total 
supply;  that  is, 

n        m 

t:  d1  =  z: 

In  this  section  a  transportation  problem  is  solved  by 
the  simplex  algorithm,  the  transportation  algorithm,  the  dis 
crete  maximum  principle,  and  the  method  of  networks. 


:  D3  =  £  °i 


An  Example.   Cities  1  and  2  are  the  supply  points  with  7 
and  3  units  of  the  product  available,  respectively.   Cities  3, 
k,  and  5  are  demand  points  with  the  demand  for  k,  k,  and  2  units 
of  the  product,  respectively.   The  shipping  costs  per  unit  are 
given  as  follows: 

To  demand  city: 
Prom  supply-    1 
ing  city:       2 

The  problem  is  to  find  the  optimal  solution  which  will  minimize 
the  total  shipping  costs. 


3 

k 

5 

5 

3 

6 

1 

k 

7 

10 


3.1.   Method  I,  Simplex  Algorithm 

The  linear  programming  problem  formulation  of  the  above 
problem  is  obtained  by  using  the  equations  (1),  (2),  and  (3) • 
The  problem  becomes  one  of  finding  x±.^   0,  the  amount  shipped 
from  source  i  to  destination  j,  i  =  1,  2  and  j  =  3,  K,    5  which 
minimizes  the  objective  function, 

z  =  5x13  +  3xxl|  +  6xl5  +  x23  +  lpc2[|_  +  7x25   [$) 
subject  to  the  supply  constraints, 


x13  +  xll|  +  xl5  =  7 
x23  +  x2l]_  +  x25   =  3 

and  the  demand  constraints, 


(6) 


x 


;13  +  x23  ~  *+ 

xll^  +  x21+  =  ^ 


(7) 


xl5  +  x25  =  2 
This  linear  programming  problem  can  be  solved  using  the 
simplex  algorithm  as  given  in  any  standard  book  on  operations 
research  (Dantzig  (2)  or  Hadley  (8)).   Refer  to  Table  1  for 
complete  solution  of  the  problem.   Note  that  each  equation  in 
(6)  and  (7)  can  be  represented  as  a  linear  combination  of  other 
equations,  owing  to  the  fact  that  supply  equals  demand.   This 
reduces  the  independent  constraint  equations  to  four;  the  fifth 
(any  one)  equation  is  redundant  and  can  be  discarded.   Thus  the 
number  of  basic  variables  in  this  type  of  problem  is  always 
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Table  1.   Simplex  Algorithm  for  the  Transportation  Problem. 


c  .= 

c .  . 

5 

3 

o 

] 

I). 

7 

M 

M 

H 

M 

J 

J- J 

°Bi 

x3 

b 

X13 

xlk 

xi5 

x23 

x2k 

x25 

Al 

A2 

A3 

*k 

b/a 

M 

Al 

7 

1 

1 

1 

0 

0 

0 

1 

0 

0 

0 

M 

Ao 

3 

0 

0 

0 

CD 

1 

1 

0 

1 

0 

0 

3/1  ->  min 

M 

A3 

Ak 

1). 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

M 

^ 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

k/1 

:r 

min 

-5 

-3 

-o 

-1 

-k 

-7 

0 

0 

0 

0 

z  = 

18M 

+2M 

+2M 

+M 

+2M 

+2M 

+M 

-:os:+Vs 

M 

Al 

7 

1 

1 

1 

0 

0 

0 

1 

0 

0 

0 

7 

1 

x23 

A3 

Aii 

3 

0 

0 

0 

1 

1 

1 

0 

1 

0 

0 

M 

1 

1 

0 

0 

0 

-1 

_2 

0 

_]_ 

1 

0 

If.  — y 

M 

fc 

0 

d) 

0 

0 

x 

0 

z 

0 

0 

1 

ir 

min 

3 

-i 

-3 

-6 

0 

-3 

-b 

0 

1 

0 

c 

z  = 

12M 

+2M 

+2M 

+M 

-M 

-2M 

-1 

M 

*1 

3 

1 

0 

1 

0 

-1 

0 

1 

0 

0 

-1 

3/1 

1 

x23 

3 

0 

0 

0 

1 

1 

1 

0 

1 

0 

0 

M 

A^ 

1 

0) 

0 

0 

0 

-1 

-1 

0 

-1 

1 

0 

i/i -> 

3 

xlk 

k 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

min 

15 

-5 

0 

-6 

0 

0 

-6 

0 

1 

0 

0 

z  = 

LM 

+2M 

+M 

-2M 

-M 

-2M 

-2M 

t 

M 

Al 

2 

0 

0 

a) 

0 

0 

1 

1 

1 

-1 

-1 

2  -*- 

1 

X23 

3 

0 

0 

0 

1 

1 

J_ 

0 

1 

0 

0 

5 

x13 

xlk 
min 

1 

1 

0 

0 

0 

-1 

-1 

0 

-1 

1 

0 

■ 

3 

1*. 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

j 

20 

0 

0 

-6 

0 

-5 

-11 

0 

-k 

t 

3 

z  = 

2M 

+M 

+  M 

-2M 

-  2K 

o 

X15 

2 

0 

0 

1 

0 

0 

1 

1 

1 

-1 

-1 

1 

x23 

3 

0 

0 

0 

1 

1 

1 

0 

1 

0 

0 

s 

x13 
xlk 

1 

1 

0 

& 

0 

-1 

-1 

c 

-1 

1 

0 

3 

k 

0 

1 

0 

0 

1 

0 

0 

0 

0 

1 

ZJ" 
°j 

min 

32 

0 

0 

0 

0 

-5 

-5 

6 

2 

-1 

-3 

z 

-M 

-M 

-M 

-M 

Note:  At  the  fourth  iteration  the  artificial  variables 
are  out  of  the  basis  and  the  optimality  criterion  for  minimum 
z  is  satisfied.   Therefore  the  solution  at  the  fourth  itera- 
tion is  the  optimal  solution  with  x13  =  1,    x-jl  =  k»  x]_$  =  2, 
and  x2-j  =  3,    and  the  minimum  cost  is  equal  to  32. 
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equal  to  (m  +  a  -  1) ,  where  m  is  the  number  of  origins  and  n  is 
the  number  of  destinations.   The  optimal  solution  obtained 
using  the  simplex  algorithm  is  the  optimal  solution  for  the 
transportation  problem  where  the  minimum  cost  is  given  by  z, 
as  shown  in  Table  1. 

It  is  observed  that  the  constraint  equations  in  any  trans- 
portation problem  have  a  simple  form.   All  the  nonzero  coef- 
ficients of  the  active  variables  x.  .  are  ones;  furthermore,  any 
given  x.  .  appears  in  two  and  only  two  of  the  constraints. 
(Refer  to  equations  (6)  and  (7).)   Because  of  these  special 
properties  of  the  constraints,  it  is  possible  to  solve  a  trans- 
portation problem  using  a  much  simpler  and  faster  method  known 
as  the  "transportation  algorithm",  which  will  be  discussed  in 
the  next  section. 

3.2.   Method  II,  The  Transportation  Algorithm 

The  example  solved  in  section  3.1*  can  be  written  in  a 
compact  matrix  form  as  shown  in  Table  2,  where  c.  .  is  the  cost 
of  shipping  one  unit  from  origin  i  to  destination  j,  the  0^  is 
the  quantity  available  at  origin  i  and  D .  is  the  requirement 
at  destination  j.   To  simplify  the  discussion,  the  algorithm  is 
explained  first  and  then  illustrated  by  solving  the  example. 

Steps  of  the  Algorithm. 

Step  1.   Assign  initial  allocations  by  one  of  the  three 
methods  given  below. 

1.  Northwest  corner  rule:   Starting  with  the  upper  left 
corner  of  Table  2,  make  the  appropriate  assignments 
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on  a  diagonal  path  to  the  lower  right  corner  depleting 
the  sources  and  satisfying  the  demand  and  having  no 
more  than  (m  +  n  -  1)  assignments. 

2.  Method  of  allocating  to  the  cell  with  minimum  ship- 
ping cost:   Determine  the  smallest  cost  for  a  cell  in 
the  tableau.   Suppose  this  occurs  for  cell  (i,  j). 
Make  the  maximum  possible  assignment  to  this  cell. 
Then  neglect  row  i  or  column  j  depending  on  whether 
source  i  is  depleted  or  demand  j  is  satisfied.   Repeat 
the  process  until  the  total  supply  is  exhausted  and 
all  the  demands  are  satisfied  and  where  no  more  than 
(m  +  n  -  1)  assignments  have  been  made. 

3.  Vogel's  method:   For  each  row,  find  the  lowest  o^j  and 

the  next  lowest  cost  c-t  in  that  row.   Compute 

( c .  .u  -  c.  .).   In  this  way,  m  numbers  are  obtained, 
i "it    ij  J 

Do  the  same  for  each  column  obtaining  n  additional 
numbers.   Choose  the  largest  of  these  (m  +  n)  numbers 
and  make  the  maximum  possible  assignment  to  the  cell 
with  the  lowest  cost  o*i    in  that  row  or  column  depend- 
ing upon  whether  this  number  resulted  from  a  row  or 
column  computation.   Cross  off  either  row  i  or  column  j, 
depending  upon  whether  the  source  i  is  depleted  or  re- 
quirement demand  j  is  satisfied,  and  repeat  the  pro- 
cess for  the  resulting  tableau.   When  the  maximum  dif- 
ference is  not  unique,  an  arbitrary  choice  can  be  made, 
and  if  a  row  and  column  constraint  are  satisfied  simul- 
taneously, cross  off  only  the  row  or  the  column,  not 
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both.  The  process  is  terminated  when  the  supply  is 
exhausted  and  the  demands  are  satisfied  and  when  no 
more  than  (m  +  n  -  1)  assignments  have  been  made. 

In  each  of  the  three  methods  the  (m  +  n  -  1)  allocations  con- 
stitute the  initial  basic  feasible  solution. 

Step  2.      Determine  a  set  of  (m  +  n)  numbers  u^,  i  =  1, 
.  .  .  ,  m,  and  v.,  j  =  1,  .  .  . ,  n  such  that  for  each  occupied 
cell  (i,  j),  Oil  =  u-  +  v..   Note  that  u.  and  v..  are  some  arbi- 
trarily  chosen  values  and  are  known  as  the  implicit  prices. 

Step  3-   Calculate  9j i   =   c- •  -  (u-  +  v  •)  for  each  unoccu- 
pied  cell  (i,  j).   9.  .  is  known  as  the  relative  cost  factor. 

Step  k.   Examine  the  values  of  9 . .  of  each  unoccupied  cell 
(i,  j);  if  all  are  positive  an  optimal  solution  has  been  found. 
If  one  or  more  are  negative,  the  solution  is  not  optimal  and  a 
change  in  the  existing  solution  is  possible  which  will  reduce 
the  cost. 

Step  5.   Find  a  path  to  a  negative  6 .  .  (generally  most 
negative)  cell  by  horizontal  and  vertical  movements  along  the 
shortest  path  of  assigned  cells. 

Step  o.      In  the  path  determined  in  Step  %,    find  the  maximum 
allocation  for  each  horizontal  move,  then  find  the  maximum 
allocation  for  each  vertical  move,  select  the  minimum  of  these 
values  and  assign  it  to  the  cell  with  the  negative  9 .  .  in 
Step  5-   Adjust  the  assignments  on  this  path  so  that  all  demands 
and  supplies  are  compatible,  then  go  to  Step  2  and  repeat  the 
process  until  an  optimal  assignment  is  found  according  to  the 
criterion  in  Step  L\..      In  solving  the  example  problem  by  this 


IS 

method,  the  initial  solution  is  found  using  the  second  method 

of  Step  1,  and  the  allocations  which  result  are:   x2->  =  3, 

x, ,  =  k,  x,  .  =  1,  and  x,,  =  2,  which  are  shown  in  Table  3.   Bote 
llj.   ^'   13  lb     ' 

that  an  occupied  or  assigned  cell  (i,  j)  is  also  known  as  the 
basic  variable  x.  .  =  1  and  the  other  cells  are  known  as  nonbasic 
variables  equal  to  zero.   Thus  cell  (2,ij.)  and  cell  (2,5)  in 


the  above  problem  correspond  to  the  nonbasic  variables  XgK  =  ^ 
and  Xpc-  =  0.   According  to  Step  2,  the  implicit  prices  u^  and 
v.  for  each  of  the  basic  variables  are  found  from  the  relation, 

cij=ui+vj  C8) 

for  the  occupied  cells. 

The  procedure  is  usually  initiated  by  setting  the  Uj_  or 
v-  with  the  most  assignments  in  the  row  or  column  equal  to  zero 
and  evaluating  the  remaining  values  by  equation  (8).   The  values 
of  Uj  and  v.  for  the  example  are  shown  in  Table  4.   The  relative 
cost  factor  6.  .  for  each  of  the  unoccupied  cells  is  evaluated 
from  equation  (9), 

eij  =  cij  -  K  +  vj>  (9) 

for  the  unoccupied  cells.   Thus 

e2ij_  =  k  -  i-k  +  3)  =5 

and  92£  =  7  -  (-li  +  6)  =5. 

Since  both  of  these  relative  cost  factors  are  nonnegative,  this 
solution  is  the  optimal  basic  feasible  solution  to  the  problem 
as  noted  in  Table  I4.. 

Summarizing,  then,  the  optimal  solution  is, 
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Table  2.   Transportation  example. 


To    3 


Destination 
k.  5      Supply 


Prom 


Origin 


Demand 


c13=5 


^r1 


4=D- 


°ll+=3 


c2ir^ 


k-=\ 


ol5=6 


7  =  o. 


3  -  0, 


025=? 


2= 


D5 


10    >z  oi  = 
i=l,2 

3=3  A,  5 


Table  3.   Initial  basic  feasible  solution. 
Destination 


To 


Prom 


Origin 


Demand 


k 


Supply 


l=x13 

5 

3 

2=xl5 

G 

7 

3=x23 

1 

^ 

7 

3 

k 

^ 

2 

10 

Table  1]..   The  optimal  solution. 
Destination 


To 


From 


Origin 


Demand 


IV 


k 


5     Supply   i 


'2V 


=5 


^5- 


=5 


u.=0 


U?=-i{. 


10 


Vi=$ 


vj,=3 


v^=6 
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Basic  variable   From  city  To  city  Ho.  of  units  transported 

x13             1          3  1 

x^           1         k  k 

x23           2        3  3 

and  the  total  minimum  cost,  z 

=  5x13  +  3x11|  +  6x15-  +  x2.  +  lpc2^  +  7x25 

=  5+12  +  12  +  3 
=  32  units  of  cost. 

3.3.   Method  III,  The  Discrete  Maximum  Principle 


The  general  algorithm  of  the  Discrete  Maximum  Principle  for 
systems  without  information  feedback  is  given  by  Fan  and  Wang 
(12)  and  is  summarized  below. 

The  process  consists  of  N  stages  connected  in  series.   The 
state  of  the  process  stream,  denoted  by  s-dimensional  vector, 

Ix  ,  x  ,  .  .  •  ,  x  1  is  transformed  at  each  stage  according 

Ll    2  sJ 


x 


to  an  r-dimensional  decision  vector  6  =  Le-i>  e2>  ■    ■    •>    erJ  » 

which  represents  the  decisions  made  at  that  stage.  The  trans- 
formation of  the  process  stream  at  the  n    stage  is  described 
by  a  set  of  performance  equations. 

n   mn  ,  n-1   n-1        n-1  '  n  „n  _n,   ,  _., 

xi  =  Ti  (x1   ,  x2   ,  ...,  x3   ;  61,  G2,  ...,  &v)       [10 1 

0 

Xj  =  ■<: 
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or  in  vector  form 

n   _n  ,  n-1    _n>    ,   0 
x  =  T   ( x    ;   9  )  and  x  =  ^ 

■where  i  =  1,  2,  .  .  .,  s;  n  =  1,    2,    .    .    .,   Kj  and  the  super- 
script n  indicates  the  stage  number. 

A   typical  optimization  problem  is  to  determine  the  sequence 
of  6n,  subject  to  the  constraints, 

,.,n  ,fin  „n         „m  <  0   n  =  1,  2,    .    .    .  ,    N     ,±1) 
^i  t9l>  92'  •  •  •'  V  *"  U>  i  =  1,  2,     .    .     .,    s     lXX; 

which  optimizes  the  objective  function  of  the  process, 

_s_      n 
S  =  5   c.  x..   c.  =  constant  (12) 

iti   1   x    1 

when  the  initial  condition  x.  =  -< .  is  given. 

The  procedure  for  finding  the  optimal  sequence  of  6'  is  to 

introduce  an  s-dimensional  adjoint  vector  z   and  a  Hamiltonian 

function  Hn  which  satisfy  the  following  relations: 

Hn  =  F-  z?  T?  (xn_1;  6n),  n  =  1,  2,    .    .    .,    N     (13) 
1=1 


z?_1  =  ~ -T  ,  n  =  1,  2,  ...,  N,  i  =  1,  2,     ...,    s  Uk) 

1      a-x?-1 


N 
z .  =  c .  ,   i  =  1,  2,    .    .    .  ,    s 

— n 
and  to  determine  the  optimal  sequence  of  control  action,  6  , 

from  the  conditions 

OEn 
H  =  extremum  or  ; =  0 . 
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Formulation  of  the  Transportation  Problem  by  the  Discrete 
Maximum  Principle.   To  formulate  the  transportation  problem  as 
a  discrete  maximum  principle  problem,  several  quantities  have 
to  be  defined.   Let  8.  =  the  quantity  of  the  commodity  sent  from 
i    source  or  origin  to  the  n   destination  or  demand  point 
where  C.  (8.)  =  the  cost  incurred  by  sending  6j.   For  a  system 

with  s  origins  and  N  destinations,  the  problem  becomes  one  of 

n 
determining  the  values  of  6 . ,  i  =  1,  2,    .    .    .,    s;  n  =  1,  2, 

.    .    . ,    N,  so  as  to  minimize  the  total  cost  of  transporting  the 

commodity, 

N    s 

csN  =  E   Z  c?  (e?) 

n=l   i=l 
subject  to  the  constraints 

(i)   6?  i   0. 


H   n 

(ii)   y__    6'.  =  0^,  the  number  of  units  of  the  commodity 
n=l 

available  at  i*-"  origin,  i  =  1,  2,    .    .    .,    s. 

i  4r      n    n 

(iii)   C-    9-  =  D  ,  the  number  of  units  of  the  commodity 

1=1 

required  by  the  n   demand  point,  n  =  1,  2,  .  .  .,  K. 

The  demand  points  are  defined  as  stages  and  the  quantity 
shipped  from  the  i    origin  to  the  first  n  stages  or  demand 
points  are  defined  as  the  state  variables  x?,  i  =  1,  2,  .  .  . , 
s  -  1.   Thus 

xi  =  0  ,  x±   =  0i 

1,  2,  .  .  .,  s-1,  n  =  1,  2,  .  .  .,  N. 


n          n-1 

n 

.    =  x.        + 

H  . 

i          i 

1 
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Note  here  that  the  number  of  units  supplied  from  the  s   -source 
to  n   destination  (stage)  can  be  obtained  by  subtracting  the 
sum  of  the  units  supplied  to  the  n    stage  from  the  first 
through  (s  -  1)    sources  from  the  total  number  of  units  re- 
quired by  the  n    stage;  that  is 

Bns   =  Dn  -  Z1  e^  (18) 

s      i=i 

A  new  state  variable,  x"  is  defined  as 

n    n-1    -  .   /.Bi      0    _ 

xs  =  xs   +  z.  c1  (e±)   ,     xs   =  o 

1=1 

i  =  1,  2,  .     .     .  ,    H 

w 
where  x'  is  equal  to  the  total  cost  of  transportation.   The 

problem  is  to  minimize  x  with  proper  choice  of  the  sequence  of 

6n,  i  =  1,  2,  .  .  .,  fa  -  1)  ;  n  =  1,  2,    .    .    .  ,  H.   The  Hamilton- 

ian  function  can  be  written  as 

n   s—   n  /  n-1    n,     n  P  n-1    .—  „n  ,„n,  '    , __, 

Hn  =  >.  z±  Ut      +  ei)  +  zs  [_xs   +  l_   ci  (ei;J   (20) 

i=l  i=l 

n  =  1,  2,  .  .  .,  N, 

and  components  of  the  adjoint  vector  are,  in  general, 


n-1    o   w      n  ,   , 

z.    =  r  =  z.  ,   i  =  1,  2,    .    .    .,    s      (21) 


with  z3  =  1,  n  =  1,  2,  .  .  .,  N. 

Since  z.  and  x.    are  considered  as  constants  at  each  step 
in  the  minimization  of  the  Hamiltonian  function  given  by  equa- 
tion (20),  it  is  convenient  to  define  the  variable  part  of  the 
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Hamiltonian  function   as 


b?  -  T   .?  e£  +  £  o;  (ej)  .  (22) 

i=l         1=1 


The  transportation  problem  previously  solved  in  sections  3-1 
and  3.2  is  presented  in  Table  $   and  will  be  solved  again  as  a 
discrete  maximum  principle  problem.   Since  in  this  example  there 
are  two  origins  and  three  destinations,  then  s  =  2  and  N  =  j. 
C1?  is  the  cost  incurred  in  shipping  one  unit  from  the  i  i  origin 
to  the  n   destination.   The  total  number  of  units  required  by 
N  destination  is  equal  to  the  total  number  of  units  supplied 
from  s  origins;  that  is, 


N         s 

X  Dn  =  £  0i 


n=l      1=1 

The  variable  part  of  the  Hamiltonian  function  using  equation 
(22)  for  the  above  example  is, 

Hv  =  zl  61  +  Cl   Gl  +  C2  92>    a  "  1.  2.  3. 

n    n    n 
Since   6?  =  D   -  8-,  ,  the  following  is  obtained. 

n.n    n    n.n    nn  „ 

Hv  =  ^zl  +  Cl  "  G2}    el  +  C2  D  '    n  =  1>    2>    i- 

Stage  1.   For  the  first  destination,  n  =  1,  the  variable 
part  of  the  Hamiltonian  equation  becomes, 

H^  =  (z\   +   C±   -   c|)  6i  +  C2  D1 

From  Table  5,  this  becomes 


Table  £•   Transportation  costs,  supplies, 
and  demands. 
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Origin 


Demand 


Destination 
2         3 


Supply 


k  -  5 

5*  -3 

4- 

.    1 

i 

°i  - 

7 

o}  =   1 

ol-k 

4- 

7 

°i  = 

3 

D1  =  k 

D2  =  k 

D3  = 

2 

i 

10 

Table 


Conditions  necessary  for  EC 
to  be  minimum. 


n\ 

n 
Minima  of  H 

occurring  at 

Sl 

n 

zi 

0 

>  -k 

1 

o  5  ej  s  k 

=  -k 

k 

<-k 

0 

>   1 

2 

0  $  Q2  %  h. 

=  1 

k 

<1 

0 

>  1 

3 

o  s  e^  $  2 

=  1 

2 

<  1 
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H*  =  ( z\  +  k)    &l   +  k   • 

Thus  z,  =  -ij.  =  C„  -  C,  is  the  boundary  point  such  that  H   is 
minimum  for  0  £  6,  $  1±. 

There  are  three  conditions  at  which  H  may  be  minimum: 

(a)  H  =  minimum  at  9-,  =  0  when  z-,  >  -ij. 

(b)  K  =  minimum  at  0  £  9.  S  l\.   when  z,  =  -1|. 

(c)  Hv  =  minimum  at  6,  =  i|.  when  z,  <  -I4.. 


Using  this  same  approach,  the  conditions  on  8,  for  minimum 
H^  are  illustrated  in  Table  6.   As  given  by  equation  (21),  the 
value  of  z£,  n  =  1,  2,  3,  is  identical.   Prom  Table  6,  the 
boundary  values  of  st^  are  determined  and  illustrated  in  Pig.  7. 

The  value  of  z!?  which  gives  all  the  solutions  satisfying 
the  constraints  from  conditions  (i),  (ii),  and  (iii)  can  be  ob- 
tained.  The  solution  which  minimizes  the  cost  is  then  selected. 
For  example,  the  solutions  6,  and  the  corresponding  values  of 
jQ   in  the  region  of  -ij.  <  z1   <  1,  n  =  1,  2,  3,  are  obtained  us- 
ing Table  6.   Knowing  the  values  of  6-,,  8,  are  obtained  from 
equation  Cl8),  they  are  as  shown  in  Table  7-   This  solution  does 
not  satisfy  the  end  point  conditions:   0-,  =  7  and  On   =  3.   In 
order  to  satisfy  the  end  point  conditions,  the  corresponding 
solution  of  z^  =  -I4.  is  found  where  0  C  9,  £  i+,  and  thus  9T  has 
to  be  equal  to  1.   By  this  process  other  solutions  are  checked 
for  feasibility  and  then  minimum  cost.   It  turns  out  8?"  =  1 


2^ 

Stage 

n=l             2 

,3 

-*— 

a 

— *- 

— >             -c — 

— 1 — 

-6  - 

-i i — 

■5    -k 

— 1 1 1 1 

-3   -2   -1   0 

1 

H 1 h- 

2         3k 

— 1 — 

5 

— 1— 
6 

-«s 

n 

i 

Pig.  7. 

Boundary  values  of  adjoint 

n 
vector  z. 

Table  7. 

0.  corresponding  to 
of  -I).  <  z"  <  1. 

value 

Destination 

n 

i 

1 

2         3 

Supply, 

0 

I 

1 

0 

k 

2 

6 
(7) 

Origin 

2 

Demand,  D 

k 

0 

0 

i   ^ 
(3) 

k 

fc 

2 

10 

Tat 

le  8. 

Optimal  feasible  solution. 

Destination 

n 
l 

1 

2        3 

Supply, 

°i 

~] 

1 
Origin 

2 

1 

k 

2 

7 

3 

0         0 

3 

10 

Demand,  Dn 

k 

|      *      |     . 
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gives  the  feasible  optimal  solution  as  shown  in  Table  8. 

"i  2 

->  -   n  n 

The  total  minimum  cost  =  >_   2_  C  •  8 .  =  5  X   1  +  3  x  Ij. 

n=l   1=1 

+6x2+1x3=  32  units,  which  is  the  same  as  the  one  ob- 
tained by  other  methods. 

3.1+.   Method  IV,  Method  of  Network 

In  this  method  a  network  representation  of  the  problem  is 
constructed  and  the  calculations  are  made  directly  on  the  net- 
work.  The  transportation  problem  considered  in  the  previous 
sections  can  be  represented  in  the  form  of  a  network  as  shown 
in  Fig.  8.   The  cost,  o. .,  of  shipping  a  unit  of  commodity  from 
origin  i  to  a  destination  j,  is  shown  on  the  respective  arc. 
The  amount  of  supply  at  each  origin,  0^,  and  the  amount  of  de- 
mand at  each  destination,  Dj,  are  noted  at  each  respective  node. 
This  type  of  problem  can  be  solved  readily  by  the  method  of 
networks,  the  steps  of  which  are  given  below. 

Steps  of  the  Method. 

Step  1.   To  find  an  initial  basic  feasible  solution  by 
forming  a  tree:   To  do  this,  each  origin  i  is  joined  to  the 
other  nodes  of  the  network  using  the  arcs  pointing  away  from 
the  origin.   Similarly,  each  destination  is  joined  to  the  other 
nodes  using  the  arcs  pointing  into  the  destination.   This  is 
continued  until  a  chain  of  arcs  from  an  origin  joins  with  a 
chain  of  arcs  to  a  destination.   At  this  point  start  joining 
other  nodes  to  this  chain  so  that  no  loops  are  formed.   When 
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all  nodes  are  joined  in  this  manner  a  tree  will  have  been  formed 
as  defined  in  section  2.   Then  assign  values  to  the  basic  vari- 
ables corresponding  to  the  arcs  of  this  tree  such  that  all  the 
supply  and  demand  conditions  are  met. 

Step  2.      Compute  the  implicit  price  %.    such  that  equation 
(23) 

"3  -  -i  =  o13  (23) 

is  true  for  arcs  (i  -t-   j)  of  the  tree  in  step  1.   Note  that 
■k^   and  1j  of  this  method  correspond  to  u^  and  v -•  of  the  trans- 
portation algorithm  in  equation  (8). 

Step  3-      Calculate  the  relative  cost  factor,  sj_-;>  using 
equation  ( 2l\.)  ,  for  each  arc  (i  — >•  j)  which  is  not  part  of  the 
tree  but  exists  in  the  original  network. 

hi   =  cij  "  (-j  -  *1>  (2^ 

Here  again  6.  .  corresponds  to  6-_.  of  the  transportation  algo- 
rithm in  equation  (9). 

Step  Ij..   Examine  the  values  of  6^-.   If  all  are  nonnegative, 
then  this  solution  is  the  optimal  basic  feasible  solution.   If 
one  or  more  6, •  are  negative,  then  the  solution  is  not  optimal; 
then  go  to  step  $. 

Step  5.   Draw  the  arc  with  a  negative  (generally  most  nega- 
tive) value  of  6 • .  and  include  that  arc  in  the  tree.   The  re- 
maining values  of  the  basic  variables  around  the  loop  thus 
formed  are  adjusted  with  a  value  assigned  to  the  new  arc.   The 
appropriate  arc  is  dropped  to  obtain  a  new  tree.   The  procedure 
is  again  repeated  from  step  2  using  the  new  basic  feasible 
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Origin 


Destination 


Pig.  8.   Network  of  transportation  problem. 
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Solution  until  the  optimal  solution  is  obtained.   The  previous 
example  as  illustrated  In  Pig.  8  is  again  solved  using  the 
method  of  networks.   In  accordance  with  step  1,  an  initial  basic 
feasible  solution  is  determined  and  is  shown  in  Pig.  9.   Note 
that  the  values  of  the  variables  assigned  on  the  arcs  of  the 
tree  in  Pig.  9  should  be  such  that  the  algebraic  sum  at  each 
node  is  zero  where  values  are  considered  to  be  positive  when 
the  arrow  points  into  the  node  and  are  considered  to  be  nega- 
tive when  the  arrow  points  away  from  the  node. 

According  to  step  2,  the  implicit  prices  are  obtained  and 
are  shown  in  Pig.  10.   In  determining  those  implicit  prices, 
the  value  it^  =  0  is  assigned  to  the  node  i  which  has  a  large 
number  of  arcs  pointing  away  from  it.   In  Pig.  9,  this  node  is 
1,  so  it-i  is  given  a  value  of  zero  and  other  values  are  obtained 
from  the  equation  (23)  which  is  3t  .  -  it.  =  c.  .,  and  are  shown 
in  Fig.  10. 

The  arcs  which  are  not  in  the  tree  but  exist  in  the  origi- 
nal network,  Pig.  8,  are  (2  — >   ij.)  and  (2  — >  $)'.      Prom  step  3, 
the  relative  cost  factor,  6.,,  for  each  of  these  arcs  is  ob- 
tained by  using  equation  (2]|),  which  is  6j  j  =  c  •  -■  -  (it*  -  It*). 

•'■J     -1- J     J     •L 

Therefore 


s2I|_  =  c21+  -  (%^  -  it2)  =  Ij.  -  (3  -  k)   =  5 

62$  =   c25  "  'K5  -  Tt2)  =  7  -  (6  -  1|_)  =5. 
According  to  step  Ij.,  since  all  o^.-'s  are  nonnegative,  this 


and 
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On  =  7 


X13  = 

1      / 

xlk  =  '+ 

>(k) 

\U 

\xi5  = 

2 

D5=2 

02=3 


Pig.  9.   Initial  basic  feasible  solution. 
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D,  =  k 


0o  =   3 


Pig.  10.   Feasible  solution  with  implicit 


prices. 
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solution  is  the  optimal  basic  feasible  solution. 

Note:  From  steps  3  and  I4.  it  is  confirmed  that  the  solution 
in  Fig.  10  is  the  optimal  basic  feasible  solution.   Thus  the 
optimal  solution  is: 

Basic  variable   From  city   To  city   Mo.  of  units  transported 
x13  1         3  1 

X^  1  k  Ij- 

x1$  l        5  2 

x23  2        3  3 

and  the  total  minimum  cost 

=  5xl   +  3xlj_+6x2   +   lx3 

=  32  units 
which  is  the  same  as  the  one  obtained  by  the  other  methods. 

3.5-   Discussion  of  Methods  of  Section  3 

The  example  transportation  problem  used  to  illustrate  the 
various  methods  is  small  in  size.   Thus  the  advantages  of  the 
various  methods  may  not  be  obvious.   In  general,  to  solve  a 
transportation  problem,  the  simplex  algorithm  takes  more  itera- 
tions than  transportation  algorithm.   The  method  which  depends 
on  the  discrete  maximum  principle  is  inefficient,  particularly 
since  the  problem  can  be  solved  with  the  simple  transportation 
algorithm.   The  problem  is  visualized  easily  using  the  method 
of  networks.   The  fastest  and  simplest  method  seems  to  be  the 
transportation  algorithm,  but  there  is  one-to-one  correspondence 
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between  this  method  and  the  method  of  networks.   In  some  class 
of  the  transportation  problems  such  as  the  transshipment  prob- 
lems, the  method  of  networks  decreases  the  number  of  quantities 
to  be  calculated  and  thus  may  actually  be  faster  than  the 
transportation  algorithm.   This  comparison  provides  the  basis 
of  the  discussion  in  section  I).. 
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k..      TRANSSHIPMENT  PROBLEM 

A  transshipment  problem  is  a  special  type  of  a  transporta- 
tion problem  where  there  are  many  possible  routes  for  shipping 
between  an  origin  and  a  destination.   Since  there  are  several 
possible  routes  from  an  origin  i  to  a  destination  J,  it  is  pos- 
sible that  there  may  be  no  direct  route  and  the  commodity  will 
have  to  be  transported  through  intermediate  points  where 
neither  supply  nor  demand  exists. 

The  various  quantities  and  terms  of  the  problem  are  de- 
fined in  the  following  manner.  Let  n  be  the  total  number  of 
nodes,  which  includes  the  origins,  the  destinations,  and  the 
intermediate  nodes.   Let 

0j_  =  net  supply  at  node  i,  i  =  1,  2,  .  .  .,  n 
D<   =  net  demand  at  node  j,  j  =  1,  2,  .  .  .,  n 

c-  •  =  c  ,-j  =  cost  of  shipping  a  unit  from  node  i 
-'-J     J 

directly  to  node  j 

°ii  =  CJj  =  ° 

x.  ■  =  quantity  shipped  from  node  i  to  node  j 

x..  =  quantity  transshipped  through  j  to  some  other  node. 

If  Cii  = «0  ,    this  implies  that  direct  shipment  between  nodes  i 

and  j  is  impossible. 

A  condition  of  the  problem  is  that  the  total  supply  equals 

the  total  demand;  that  is, 


i=i    x     j=i    ° 


J4 


Thus  the  problem  is  one  of  determining  the  quantity  of  the 
commodity  supplied  from  origin  i  to  destination  j,  x^.^0, 
which  minimizes  the  objective  function, 


z 


n    n 
V 


i_   Zl  C.J.-  Xjj,   where  i  ^  j  (26) 


i=l   j=i 
subject  to  the  supply  constraints, 

t    xii  "  xii  =  °i>   i  =  1,  2,  .  .  .,  n      (27) 
j=l    J 

i  f'  j 

which  states  that  the  total  amount  shipped  out  of  node  i  less 

amount  through  node  i  equals  net  supply  from  node  i  and  the 

demand  constraints, 

jj^ij  -  Xjj  =Dj>   J  =  1,  2,  .  .  .,  n      (26) 

i  71   J 
which  states  that  the  total  araount  shipped  in  node  j  less 
amount  through  node  j  equals  net  demand  at  node  j. 

It  is  noted  that  there  are  (2n  -  1)  basic  variables  x. ,. 
An  example  problem  consisting  of  eight  nodes,  two  source 
points,  three  demand  points,  and  three  intermediate  points,  is 
given  below.   It  will  be  set  up  as  a  linear  programming  problem, 
and  will  be  solved  by  both  the  transshipment  algorithm  and  the 
method  of  networks. 

An  Example.   There  are  \$   lots  of  canned  tuna  in  Seattle 
and  30  in  Los  Angeles.   Of  these,  2£  are  to  be  sent  to  Boston, 
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35   to  New  York,  and  15  to  Miami.   The  shipping  cost3  per 
between  the  corresponding  cities  are  illustrated  below: 


lot 


Denver 

Kansas 

City 

G  h  i  o  t: 

go 

B 

oston 

New 

York 

Miami 

Seattle 

11 

•::- 

20 

Los  Angeles 

10 

15 

20 

* 

* 

Denver 

0 

5 

10 

25 

* 

Kansas  City 

5 

0 

10 

15 

Chicago 

10 

10 

0 

13 

12 

11 

New  York 

12 

5 

0 

'"'No 

direct 

r> 

oute  available. 

The  problem  is  one  of  finding  the  optimal  shipping  routes 
which  will  minimize  the  total  shipping  costs. 

Before  discussing  the  different  methods  for  solving  this 
problem  it  is  helpful  to  visualize  the  problem  as  the  network 
problem  shown  in  Pig.  11. 

I4..I.   Method  I,  Simplex  Algorithm 

As  in  the  case  of  the  transportation  problem,  using  equa- 
tions (2b),    (27),  and  (28),  the  transshipment  problem  can  be 
formulated  as  a  linear  programming  problem  and  solved  by  the 
simplex  algorithm.   The  problem  as  formally  stated  becomes  one 
of  determining  the  x.  .  J-  0  which  minimize  the  objective  function, 


i=l   j=l 

i  ¥  i 


c .  .  x.  .,  where  i  =  1,  2, 
J  =  1,  2, 


(29) 
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I.N. 


I.N. 


Orn.  :  Origin 

Dest.:  Destination 

I.N,  :  Intermediate  node 


1.  Seattle 

2.  Boston 

3.  Denver 
I4..    Chicago 


5.  New  York 

6.  Kansas  City 

7.  Los  Angeles 

8.  Miami 


Pig.  11.   Network  of  transshipment  problem 
with  c-  ,,    0, ,    and  D~-. 
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subject  to  the  supply  constraints, 

-xu  +  x12  +  x13  +  xlh   +  xl5   +  xl6  +  x1?  +  xl8  =  kS 

x21  "  x22  +  x23  +  x2li  +  x25  +  x26  +    x27  +  x28  =  ^ 
(30) 

x71  +  x72  +  x73  +  x7l+.  +  x75  +  x?6  "  x77  +  x78  =  3° 

x8i  +  x82  +  x83  +  x81^  +  x85  +  x86  +  x87  "  x88  =  ° 
and  the  demand  constraints, 

"xil  +  x2i  +  x3l  +  xl|l  +  x5l  +  x6l  +  x7l   +  x8i  =  ° 

x12  "  x22  +  x32  +  \2   +  x52  +  x62  +  x72  +  x82  =  2^ 

(3D 

xi5  +  x25  +  x35  +  xk5  '   x55  +  x65  +  x75  +  x85  =  35 

xl8  +  x28  +  x38  +  xk8  +  x58  +  x68  +  x78  "  x99  =  l5 

The  total  number  of  dependent  equations  in  (30)  and  (31) 
are  l6,  and  eliminating  any  one  equation  leads  to  a  set  of  15 
independent  constraint  equations.   Thus  the  linear  programming 
problem  reduces  to  a  problem  with  1$   equations  and  79  variables 
(61j.  active  variables  +  15"  artificial  variables). 

The  problem  given  in  this  form  (equations  (29),  (30),  and 
(3D)  is  easily  solved  using  simplex  algorithm  as  was  shown 
for  the  transportation  problem  in  section  3-1  (Dantzig  (2)  or 
Hadley  (8)).   This  problem  was  solved  on  an  IBM  1620  computer 
and  the  optimal  solution  was  obtained  after  19  iterations. 
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1)..2.   Method  II,  The  Transshipment  Algorithm 

This  method  is  precisely  the  transportation  algorithm 
explained  in  section  3.2,  but  is  known  as  the  transshipment 
algorithm  when  used  to  solve  transshipment  problems.   The  only 
difference  is  that  in  the  transshipment  algorithm  an  initial 
basic  feasible  solution  is  found  by  forming  a  network  tree,  as 
noted  in  step  1  of  section  3  •  h>    in  lieu  of  the  method  suggested 
in  step  1  of  section  3-2. 

The  steps  of  the  algorithm  can  be  summarized  as  follows: 

Step  1.   Find  an  initial  basic  feasible  solution  forming 
a  tree  (refer  to  step  1,  section  3-k-)>    an<3  form  a  tableau. 

Step  2.      Determine  a  set  of  2n  arbitrary  numbers,  Uj  and 
v.,  such  that  for  each  occupied  cell  (i,  j),  cn-  ■  =  u^  +  v ,. 
Note  that  these  arbitrary  numbers  u-,  v-  are  the  same  implicit 
prices  as  before. 

Step  3-   Calculate  9----,  which  is  known  as  the  relative 
cost  factor,  for  each  of  the  unoccupied  cells  (i,  j)  such  that 

9ij  =  cij  "  <ui  +  -j'- 

Step  Ij..   Examine  the  values  of  G.  . ;  if  all  are  nonnegative, 

then  this  solution  is  optimal.   If  one  or  more  G.  •  are  negative, 

then  this  solution  is  not  optimal  and  a  change  in  the  existing 

solution  is  possible  which  will  reduce  the  cost. 

Step  5-  Find  a  path  to  the  cell  which  has  a  negative  9.  . 
(generally  most  negative)  by  horizontal  and  vertical  movements 
along  the  shortest  path  of  occupied  cells. 

Step  6.   In  the  path  determined  at  step  5  find  the  maximum 
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assignment  for  each  horizontal  move,  then  the  maximum  assign- 
ment for  each  vertical  move  and  select  the  minimum  of  these 
values,  and  assign  it  to  the  cell  with  the  negative  e.  .  in 
step  5. 

Step  7.   Adjust  the  assignments  on  this  path  so  that  all 
demands  and  supplies  are  compatible. 

Step  8.   Go  to  step  2  and  repeat  the  procedure  until  an 
optimal  assignment  is  found  according  to  the  criterion  in  step  L. 

In  solving  the  example  problem  by  this  method,  the  initial 
basic  feasible  solution  is  found  according  to  step  1,  as  shown 
in  Pig.  12.   Note  that  the  tree  evaluates  (n  -  1)  basic  vari- 
ables out  of  a  possible  (2n  -  1)  basic  variables,  where  n  is 
the  total  number  of  nodes  present  in  the  problem.   Thus  in  the 
above  example  where  n  equals  8,  seven  basic  variables  are  ob- 
tained from  the  tree  in  Fig.  12  and  the  other  eight  basic  vari- 
ables are  the  transshipment  variables  themselves,  that  is,  x.,, 
Xoo,  •  •  •>  xfl8-   Prom  this  feasible  solution  the  elements  of 
the  transshipment  algorithm  are  examined  and  are  shown  in 
Table  9. 

The  implicit  prices  u^,  v.-  and  the  relative  cost  factors 
6j  •,  for  the  unoccupied  cells,  are  determined  and  are  illus- 
trated in  their  respective  positions  in  each  cell  of  Table  9. 
In  the  transshipment  problems  note  that  the  transshipment  basic 
variables,  which  are  x-,-.,  X22,  ■    ■    ■,    X8Q»  carL  ^e  negative  since 
they  do  not  contribute  to  the  objective  function  because  of 
their  zero  cost  coefficients,  c, .  =  c„o>  •    ■    •>    -   OrS  =  ®'      ^^us 
it  is  seen  from  Table  9  that  x-,-,,  Xll,  and  x^  are  negative 


^0 


15 


Fig.  12.   Initial  basic  feasible  solution 
( step  1) . 


kl 


Table   9.      Transshipment   tableau,    initial 
basic   feasible    solution. 


DESTINATION 


Dera 


k2 


quantities.   It  is  observed  that  u<  =  -Vj  when  c< ,  =  0,  where 
i  =  j  for  each  of  the  transshipment  variables. 

Observing  in  the  above  example  that  two  of  the  relative 
cost  factors  are  negative,  they  are  9-7^  and  &*>&,    thus  according 
to  step  ij.  this  solution  is  not  the  optimal  solution.   Since 
both  S-jz    and  @j^   are  equal  to  -1,  either  variable  may  be 
selected  to  enter  the  basis;  in  this  instance  Xy^  is  selected. 
Continuing  at  step  5,    a  path  to  the  cell  (7,  6)  is  drawn  by 
horizontal  and  vertical  movements  along  the  shortest  path  of 
occupied  cells,  as  shown  in  Table  9. 

According  to  step  6,  the  value  of  Xy^  to  be  assigned  is 
determined  to  be  15,  and  the  variable  Xo£  is  dropped  to  maintain 
feasibility  of  supply  and  demand.   This,  then,  is  a  new  basic 
feasible  solution  which  is  checked  for  optimality  by  repeating 
the  procedure  from  step  2.   The  second  iteration  for  the  ex- 
ample is  illustrated  in  Table  10.   The  basic  solution  is  not  yet 
optimal  as  it  is  seen  that  the  relative  cost  factor  for  Xj~    is 
negative.   After  drawing  the  stepping  stone  path  it  is  observed 
that  the  only  change  possible  is  the  interchanging  of  the  arc 
(1  — s-  3)  and  (7  — >   3).   Thus  the  new  basic  variable  to  enter  is 
Xyn  which  will  have  a  value  of  zero  since  the  variable  to  leave 
the  basis  is  x,  ,  =  0,  as  shown  in  Table  10.   The  only  change 
in  implicit  prices  is  for  u,  and  v,  which  are  -10  and  10,  re- 
spectively.  The  resulting  basis  thus  formed  is  the  optimal 
solution  since  all  of  the  G • •  values  are  nonnegative.   Summariz- 
ing, then,  the  optimal  solution  is, 


Table   10.      Transshipment   tableau,    second 
iteration. 


DESTINATION 
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Dem.    D 


kk 


Basic  variable 
(those  contribute 
to  the  objective 

function) 

X7k 
x76 

*k5 

x68 


Prom  city 


To  c i ty 


Seattle  Chicago 

Los  Angeles  Denver 

Los  Angeles  Chicago 

Los  Angeles  Kansas  City 

Chicago  Boston 

Chicago  New  York 

Kansas  City  Miami 


No.  of  lots 
transported 


kS 

o 

IS 

25 

35 
15 


and  the  minimum  shipping  cost  =  z 

=  k5   x  20  +  15  x  20  +  15  x  15  +  25  x  13  +  35  x  12 
+  15  x  15  =  2395  units. 

I4..3.   Method  III,  Method  of  Networks 


In  this  section  the  transshipment  problem  is  solved  employ- 
ing the  method  of  networks  explained  in  section  3 . it-.   This 
method,  of  course,  has  a  one-to-one  correspondence  with  the 
transshipment  algorithm  but  in  lieu  of  using  a  matrix  format 
the  quantities  are  directly  calculated  on  the  network  graph 
itself. 

The  steps  of  the  method  can  be  summarized  as  follows: 

Step  1.   Find  an  initial  basic  feasible  solution  forming 
a  tree  (refer  to  step  1,  section  3. 1+)  . 

Step  2.   Compute  the  implicit  price  %■    such  that  Jtj  -  itj 
=  c-  ■  is  true  for  each  arc  (i— »  j)  of  the  tree  in  step  1. 
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selected).  The  values  of  the  basic  varia 
are  adjusted  with  a  value  asaigm  -  ;he 
basis.   The  appropriate  arc  is  dropped  and   new  tree  is  formed 

Ste£  6.   Using  the  new  basic  feasible  solution,  the  pro- 
cedure is  repeated  from  s^ep  2  until  zhe   optimal  solution  is 
obtained.   The  transshipment  problem  which  was  treated  in  sec- 
tions k.l   and  k-2   is  again  illustrated,  below  using  the  method 
of  networks.   The  example  consists  of  eight  nodes  which  in- 
cludes two  origins,  three  destinations,  and  three  intermediate 
points.   The  cost  of  shipping  the  demand  requirements  and  the 
quantities  available  are  as  shown  in  Pig.  11. 

An  initial  basic  feasible  solution  including  the  implicit 
prices  and  the  relative  cost  factors  from  steps  1,  2,    and  3  is 
shown  in  Pig.  13.   It  is  noted  in  Pig.  13  that  the  optimality 
criterion  is  not  met  with  the  initial  basic  feasible  solution 
because  S?3  =  -l.   This  indicates  that  it  is  profitable  to 
transport  through  arc  (7->3).   Nevertheless,  when  arc  (?-.  3) 
is  introduced,  the  subnetwork  does  not  construct  a  tree  since 
it  contains  a  loop  shown  in  Pig.  lk.      Since  the  values  of 
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Pig.  13.   Initial  basic  feasible  solution  with 
implicit  prices,  3Cj,  and  relative 
cost  factors,  5_.  ■ . 
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variables  constituting  the  loop  are  affected  by  shipping  through 
arc  (7-*-  3) ,    we  need  only  consider  these.   From  the  loop  it  is 
seen  that  the  maximum  allowable  value  of  new  variable  30?,  is 
zero,  since  the  minimum  value  of  the  arc  with  negative  sign  is 
zero.   To  avoid  the  problem  of  the  loop  and  to  allow  flow 
through  (7->3),  the  arc  (1— «-3)  is  to  be  dropped.   This  then 
forms  a  new  basic  feasible  solution,  as  is  shown  in  Pig.  15>. 
The  optimality  criterion  is  met  with  this  solution  since  all 
the  %■' s   remain  the  same  and  all  the  other  &-^'s  are  nonnega- 
tive.   Therefore  this  solution  to  example  in  Pig.  11  is  the 
optimal  basic  feasible  solution,  and  is  the  same  as  the  one 
obtained  in  section  l\..2. 

lj..I|_.   Discussion  of  Methods  of  Section  I4. 

In  solving  the  transshipment  problem  by  the  various 
methods,  it  is  seen  from  this  example  and  which  is  generally 
true  that  the  simplex  algorithm  requires  more  variables  and  mora 
iterations  than  the  transshipment  algorithm.   To  solve  a  prob- 
lem manually  the  method  of  network  has  several  advantages  over 
the  simplex  algorithm  and  the  transshipment  algorithm.   Pirst 
the  number  of  quantities  to  be  calculated  are  always  less  using 
the  method  of  network.   Secondly,  it  is  easy  to  visualize  the 
problem  through  the  network.   Specifically,  comparing  the 
methods  applied  to  the  example  considered  in  this  section,  it 
is  seen  that  the  linear  programming  problem,  where  the  simplex 
algorithm  can  be  applied,  consisted  of  15  equations  and  79 
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Oi  =  K>- 


*73 


•*(£>■  -^.-15  \] 


o?  =  30 


Pig.  llj..   To  find  the  value  of  new  basic 
variable,  x73- 


On  =  l4&" 


Dn;  =  35 


Dfl  =  15 


Fig.  15.   The  optimal  basic  feasible  solution. 
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variables  (6I4.  active  variables  plus  1$   artificial  variables) 
when  the  total  number  of  nodes  is  8.   It  took  19  iterations 
to  obtain  the  optimal  solution  to  this  problem,  working  on  an 
IBM  1620  computer.   Employing  the  transshipment  algorithm,  the 
number  of  quantities  to  be  evaluated  at  each  iteration  is  36, 
and  the  optimal  solution  required  three  iterations,  whereas 
working  the  same  example  with  the  method  of  networks  there  are 
only  II4.   quantities  to  be  calculated  at  each  iteration  and  the 
optimal  solution  required  one  iteration.   Thus  it  appears  as 
the  size  of  problem  increases,  the  method  of  network  signifi- 
cantly reduces  the  number  of  calculations  in  obtaining  the 
optimal  solution. 
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5.   MAXIMAL  PLOWS  IN  NETWORKS 

The  network  problems  previously  considered  had  infinite 
arc  capacities  where  the  objective  was  to  minimize  cost  or  maxi- 
mize profit.   The  latter  is  accomplished  by  considering  the 
negative  values  of  the  profit  coefficient  at  each  arc.   The 
problems  considered  in  this  section  have  finite  arc  capacities 
and  consist  of  a  single  source  and  sink  connected  by  several 
intermediate  nodes,  as  is  illustrated  in  Pig.  16.   In  problems 
of  this  type  the  source  is  designated  node  1  and  the  sink  as 
node  n.   Thus  in  Pig.  16,  node  1  is  the  source  and  node  $   is 
the  sink.   The  flow  may  represent  such  things  as  fluids,  elec- 
tricity, and  automobiles  with  physical  units  per-  unit  time. 
The  flow  capacities,  d- •  ,  in  each  direction  are  designated  along 
each  arc  near  the  node  at  which  the  flow  originates.   The  flow 
is  assumed  to  be  steady  state,  that  is,  the  flow  does  not 
change  with  time. 

It  is  convenient  at  this  point  to  define  some  necessary 
terms  and  quantities. 

Let  x^ .-  >,  0'  denote  the  quantity  of  flow  from  node  i  to 
node  j  where  the  following  capacity  constraints  hold: 

0  &  x.  •  £  dx. 

1J     ij  (32) 

0  S  Xji  S  dj± 
and  the  conservation  of  flow  equations  hold  at  each  node  k: 


2-  xik  "  *-    xki  =  °  (33) 


with  k  =  2,  3,  .  .  .,  n  -  1. 
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Equation  (33)    states  that  the  sura  of  flows  into  node  k  is 
balanced  by  the  sum  of  flows  from  k  and  where  node  k  is  not 
the  source  or  the  sink. 

Prom  Pig.  lo  and  equation  (32),  the  following  relations 
are  true : 

0  Gx12  S  5,  0  -5  Xgl  ^  °> 

0  S  ^23^  "O  ,  0  ^  Xm  S  2, 

and  if  there  is  a  flow  of  2  through  arc  ( 2  — =- 1|_)  and  flow  of  3 
through  arc  (3-*k)*    then  the  flow  through  arc  (1+— *■  5)  must  be 
2+3=5-   Similarly,  if  f  is  the  flow  into  the  source  from 
outside  the  network,  then  by  definition  the  net  flow  at  the 
source  is, 

f  +  E  xu  -  E  Xji  =  0  (3k) 

i         J 

and  consequently  at  the  sink 


h*   xni  -  f  =  0  (35) 

J 


xin 


The  equations  (33),  (3k)  >    and  (3^>)    correspond  to  the  flow 
in  the  arcs  at  the  nodes  in  the  network. 

Maximal  Plow  Problem.   There  is  a  class  of  flow  problems 
which  are  governed  by  the  above  relationships  where  the  objec- 
tive is  to  obtain  the  maximal  flow  through  the  network.   In  par- 
ticular, the  maximal  flow  problem  is  one  of  determining  Xj_  .•  'j.   0, 
the  flow  through  arc  (i-v-j),  so  that  the  input  f  at  the  source 
is  a  maximum  and  where  equations  (32),  (33),    (3k)  >    an°  (35) 
govern  the  system.   A  feasible  flow  or  flow  path  in  a  network 
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2   3 


Source 


Sink 


Pig.  16.   A  network  with  directed  arc  capacities. 


Source'  2 


Pig.  17.   A  maximal  flow  example. 
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is  said  to  be  maximal  if  the  flow  from  source  to  sink  is  finite, 
and  where  no  other  feasible  flow  or  flow  path  yields  a  larger 
flow  from  source  to  sink.   There  are  several  methods  for  solv- 
ing the  maximal  flow  problem  and  three  are  discussed  in  the 
following  sections. 


5.1.   Method  I,    Chains  of  Positive  Arc 
Capacities  Joining  Source  to  Sink 


The  first  method  of  solving  these  problems  is  one  which 
utilizes  the  chains  of  positive  arc  capacities.   The  method 
consists  of  finding  the  chains  connecting  the  source  to  the 
sink  which  have  positive  arc  capacities.   This  method  is  direct 
application  of  the  following  theorem: 

Theorem  2.   If  a  maximal  flow  exists  and  the  capacities  are 
integers,  the  method  will  construct  only  a  finite  number  of 
positive  chain  flows,  whose  algebraic  sum  is  the  maximal  flovj 
(2).   The  steps  of  this  method  are  summarized  as  follows. 

Step  1.   Join,  the  arcs  of  positive  capacity  until  a  chain 
of  arcs  connects  the  source  to  sink. 

Step  2.  Determine  the  smallest  capacity  f  of  an  arc  on 
the  chain  found  in  step  1,  then  f  is  the  flow  in  this  chain, 
contributing  to  the  final  flow  from  source  to  sink. 

Step  3.      Decrease  the  capacity  of  each  directed  arc  on  the 
chain  in  step  1  by  amount  f  found  in  step  2,  and  increase  -Che 
same  amount  in  the  opposite  direction. 

Step  !>..   The  process  is  repeated  for  the  network  formed  in 
step  3  by  starting  again  at  step  1. 


5k 


Step  £.   This  procedure  is  repeated  until  no  more  chains 
of  positive  arc  capacity  joining  source  to  sink  are  possible. 
The  sum  of  the  flows  in  each  chain  obtained  in  step  2  is  the 
maximal  flow  in  the  network.   This  method  is  illustrated  by  trie 
following  example. 

An  Example.   The  problem  considered  is  illustrated  as  a 
network  in  Pig.  17.   The  capacities  of  the  arcs  in  each  direc- 
tion are  indicated  on  the  network.   The  problem  is  one  of  find- 
ing the  maximal  flow  through  the  network. 

The  method  of  solution  is  begun  in  accordance  with  step  1 
by  obtaining  a  chain  with  positive  arc  capacity  joining  the 
source  to  the  sink.   This  is  facilitated  by  forming  a  tree  of 
arcs  with  positive  capacities  and  selecting  the  appropriate 
chain.   One  possible  tree  is  shown  in  Pig.  18,  which  leads  to 
the  chain  (l->2),  ( 2  -*  I4.)  ,  (lj.-»-6),  with  capacities  (3,  ii,  6). 
As  directed  by  step  2,  the  minimum  capacity  among  the  arcs  of 
this  chain  is  determined  to  be  3  for  arc  (1-^2).   Thus  the  flow 
initiated  along  this  chain  is  f,  =  3,  which  is  shown  at  the 
middle  of  the  arcs;  the  arrows  indicate  the  direction  of  flow 
along  the  chain  (refer  to  Fig.  18) .   In  step  3,  the  capacity  of 
each  arc  is  adjusted  by  subtracting  f ,  =  3  from  the  capacity  at 
the  base  of  the  arrow  and  adding  it  to  the  capacity  at  the  point 
of  the  arrow.   This  result  is  expressed  in  Fig.  19.   The  pro- 
cedure is  repeated  again  for  the  new  network  starting  at  step  1. 
The  chain  of  arcs  with  positive  capacity  in  the  network  of 
Fig.  19  is  (1— >3),  (3-*W,  (!).->  6),  with  capacities  (2,  2,    3). 
Thus  the  flow  f2  =  2  is  permissible  along  this  chain.   Continuing 


55 


Source 


Pig.  18.   Initial  tree  with  positive 
chain  flow. 


Pig.  19.   Adjusted  arc  capacities  and 
s  e  c  ond  cha  in . 
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the  capacities  are  adjusted  again  which,  results  in  the  final 
network  as  shown  in  Pig.  20.   The  procedure  is  stopped  because 
no  additional  chains  of  positive  arc  capacity  connecting  1  and 
6  exist.   Finally,  the  total  maximal  flow  in  the  network  is  the 
algebraic  sum  of  the  chain  flows  as  obtained  from  Pigs.  l8  and 
19,  and  illustrated  in  Pig.  21.   Thus  the  maximal  flow 

£  *  tx  +  f2  =  3  +  2  =  5. 

The  above  method  and  the  methods  which  follow  in  this  sec- 
tion are  based  on  the  max-flow  min-cut  theorem  of  Ford  and  Ful- 
kerson  (3,  7,  8).  All  of  these  methods  lead  us  to  the  value  of 
minimum  cut  which  is  equal  to  the  maximal  flow.  Before  stating 
the  theorem  the  following  concepts  are  defined. 

Definition  of  a  Cut.   A  cut  is  any  set  of  directed  arcs 
containing  at  least  one  arc  from  every  chain  of  positive 
capacity  joining  the  source  to  the  sink. 

Definition  of  the  Cut  Value.   The  cut  value  is  the  sum  of 
the  capacities  of  the  arcs  of  the  cut. 

With  these  definitions  the  theorem  is  stated  as  follows. 

Theorem  3,  Max-flow  Min-cut  Theorem.   For  any  network  the 
maximal  flow  value  from  the  source  to  the  sink  is  equal  to  the 
minimal  cut  capacity  of  all  cuts  separating  the  source  from  the 
sink;  that  is,  max-flow  value  equals  the  min-cut  value. 

In  the  above  example,  see  Pig.  21,  the  set  of  directed 
arcs,  (l->2),  (1-^3),  with  capacities  (3,  2),  constitute  a 
cut.   The  value  in  this  case  is  3  +  2  =  £.   This  is  the  only  cut 
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Fig.  20.   Pinal  network,  no  chain  of 
positive  capacity. 


Cut  value  =  5 


Sou*rc 


Pig.  21.   The  maximal  flow 
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in  this  particular  case;  therefore  the  minimum  cut  value  is  > 
which  is  also  maximal  flow  value  as  seen  above. 

5.2.   Method  II,  The  Labeling  Technique 

The  labeling  technique  which  is  used  with  the  maximal  flow 
problems  is  a  more  systematic  way  of  determining  whether  a  flow 
in  the  network  exists  or  not.   Let  the  excess  capacity  through 
arc  ( i  — - -j)  be  defined  by 

Sij  =  dij  "  *ij  >  0  (36) 

where  <i±.    is  the  capacity  of  arc  (i— *-j)  and  x±1  is  the  net  flow 
through  arc  ( i  ---  j)  from  node  i  to  node  j.   Similarly, 

gji  =  dji  +  Xij'  Sin0e  Xji  =  "xij       (37) 
where  x.  .  may  have  a  value  at  the  subsequent  steps  of  the  pro- 
cedure but  it  is  assumed  to  be  zero  initially.   Having  defined 
these  quantities,  the  discussion  of  the  method  will  now  proceed. 
The  steps  of  the  method  are  as  listed  below: 

steP  I-   Start  by  considering  all  the  nodes  which  are  joined 
to  the  source  by  arcs  of  positive  excess  capacity.   On  the  net- 
work, label  a  node  j  with  two  numbers  (8-,  r.),  where  6.  =  g-,  ., 
and  rj  =  1-   (Note  that  the  value  of  6,  is  the  excess  capacity 
from  the  source  to  node  j,  and  the  value  of  r-  indicates  the 
node  from  where  one  started  and  proceeded  to  label  j;  in  the 
above  case  this  node  is  the  source,  therefore  r*  =  1.)   If  the 
sink  n  is  labeled,  go  to  step  I4.;  if  not,  go  to  step  2. 


Step 


9i 


2.   When  the  sink  is  not  labeled  in  step  1,  let  j  be 


the  smallest  index  of  the  node  from  the  set  of  nodes  labeled 
at  step  1.   Determine  if  there  are  any  unlabeled  nodes  which 
are  joined  to  j  by  arcs  of  positive  excess  capacity.   If  there 
are  no  such  nodes,  move  to  the  next  lowest  index  j  and  repeat 
the  process.   If  some  unlabeled  nodes  can  be  reached,  using  the 
index  k  for  these  unlabeled  nodes,  label  each  as  follows: 

5k  =  min(gjk,  6j),   rk  =  j  (38) 

The  label  6k  on  node  k  gives  the  minimum  excess  capacity  of  the 
two  arcs  which  form  the  path  from  the  source  to  j  and  from  j 

he  label  r,  indicates  the  node  from  which  one  proceeded 


to  k.   T 


lk 


to  label  k.   After  having  labeled  all  possible  nodes  k  with  the 
lowest  index  j,  if  sink  is  labeled  go  to  step  tj.j  if  not,  go  to 

step  3< 

Step  3.   Select  the  node  with  lowest  index  from  the  set  of 
nodes  denoted  by  k  in  step  2;  look  for  the  unlabeled  nodes  that 
are  joined  to  this  node  by  arcs  of  positive  excess  capacity. 
Let  these  nodes  be  denoted  by  the  index  q.   Label  the  nodes  q 
in  the  similar  way  as  before,  that  is, 

5q  =  mln(Skq>  5k>>   rq  =  k  (39) 

This  process  is  repeated  until,  in  a  finite  number  of  steps, 
one  reaches  one  of  the  two  following  states: 

(1)  Ho  additional  nodes  can  be  labeled,  and  the 
sink  is  not  labeled. 

(2)  The  sink  is  labeled. 


60 


If  state  (1)  is  reached,  the  existing  flow  is  maximal.   If 
state  (2)  is  reached,  the  existing  flow  can  be  increased;  go 
to  step  1+. 

Step  Ij..   If  5   is  the  label  on  the  sink,  then  determine 
the  chain  with  this  positive  flow  which  has  been  followed  in 
moving  from  the  source  to  the  sink.   Work  backwards  and  trace 
the  path,  noting  the  second  label  on  the  nodes  which  indicates 
the  preceding  node. 

Step  5-   Let  g. .  be  the  excess  capacities  of  the  arcs  in 
the  chain  which  lead  to  the  labeling  of  the  sink,  in  the  direc- 
tion in  which  one  moved  in  going  from  source  to  sink.   Then, 
consider  the  new  network  whose  excess  capacities  are, 


D1  .'  .'  =   O"  .  .   -   O 

<=ij   &ij   un 

for  each  arc  (i'  -*-  j '  )  in  the  chain  in  step  I4. 
Si'j'  =  gij>      gji  =  gj'i 


(ko) 


(14-D 


for  arcs  not  in  the  above  chain. 

Step  6.   Repeat  the  entire  labeling  process  considering 
the  new  network  obtained  in. step  5- 

Step  1.      Add  all  6   obtained  in  step  ij.  to  compute  the 
maximal  flow. 

After  a  finite  number  of  steps,  one  must  reach  state  (1) 
of  step  3,  since  6   is  always  strictly  positive,  which  assumes 
that  there  does  exist  a  maximal  flow. 

The  example  network  problem  considered  in  section  5-1  will 
now  be  solved  by  the  labeling  technique. 
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Assuming  initial  flow  x^  =  Xjj.  =  0,  from  equations  (36) 
and  (37),  the  excess  capacities  for  the  arcs  of  the  network 
are  illustrated  in  Fig.  23.   Since  x^  =  Xjj  =  0,  then  &± 
According  to  step  1,  the  nodes  joined  to  the  source  by  arcs  of 
positive  capacity  are  2  and  3,  as  g±2   =  3  and  g13  =  2  •   Tliese 
nodes  are  labeled  with  (3,  D  and  (2,  1) ,  respectively;  the 
first  number  indicates  6j  and  the  second  number  r..      In  this 


j  "  Gi- 


case  r?  =  r,  =  1  since  they  are  reached  from  the  source.   Ac- 
cording to  step  2,  the  lowest  index  node  from  the  nodes  labeled 
above,  that  is,  between  2  and  3,  is  node  2.   The  nodes  which 
are  not  labeled  and  are  joined  to  node  2  by  arcs  of  positive 
excess  capacity  are  I4.  and  $.      Node  3  is  ignored  as  it  has  been 
already  labeled.   Nodes  k  and  5  are  labeled;  for  example,  for 


k,  6,  =  min(62,  g2  ^)  =  min(3,  k)  =  3  and  r^ 


2 


arcs  of  positive  excess  capacity.   Therefore  according 


the  node 

since  node  k  is  reached  from  node  2  (refer  to  Pig.  23). 

sink  is  still  unlabeled;  hence  the  second  lowest  index  node 

among  those  labeled  at  step  1  is  considered.   This  is  node  3, 

but  there  are  no  unlabeled  nodes  which  are  joined  to  node  3  by 

0 

step  3,  the  lowest  index  node  from  those  labeled  at  second  step 
is  considered,  that  is,  the  node  k.   The  only  unlabeled  node 
connected  to  the  node  k  by  an  arc  of  positive  capacity  is  node 
6  which  has  a  label  (3,  k) •      Tnus  tne  sink  is  labeled;  there- 
fore the  flow  in  the  network  can  be  increased  from  zero  by 
5   ,  =  3  (refer  to  Fig.  23) .   The  chain  is  traced  working  back- 
wards as  instructed  by  step  'k:  the  sink,  node  6,  is  reached 
from  node  k,  node  k  is  reached  from  node  2,  and  node  2  is 
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Source  2 


Pig.  22.   Maximal  flow  problem  with. 
arc  capacities. 


(3M 


Fig.  23.   Maximal  flow  problem  with  excess 
capacities  and  nodes  labeled. 
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reached  from  node  1,  which  is  the  source.   Therefore  the  chain 
is  (1— *2)j  ( 2  — >  1+)  ,  (If-*- 6).   Now  utilizing  equations  (I4.O)  and 
( J4.I)  ,  capacities  of  the  arcs  are  adjusted  and  the  new  network 
is  formed  as  shown  in  Fig.  21)..   The  whole  procedure  from  step  1 
is  repeated  again.   Labels  on  the  nodes  are  also  shown  in 
Pig.  21)..   The  flow  in  the  network  is  again  ,  increased  by 
6   /  =  2  through  the  chain  (l->3),  ( 3  ~"  W  ,  (!).->  6) ,  as  illus- 
trated in  Fig.  21)..   At  this  point  the  total  flow  is  equal  to 
3+2=5-   Again  the  necessary  adjustments  in  the  capacities 
are  made  which  results  in  new  network  as  shown  in  Fig.  25. 
It  is  observed  that  no  node  can  be  labeled  from  the  source, 
node  1,  and  thus  the  sink,  node  6,  is  unlabeled.   This  indi- 
cates that  state  (1)  of  step  3  is  reached,  and  consequently  the 
existing  flow  is  maximal.   Thus  the  maximal  flow  for  this  ex- 
ample is  5,  which  is  the  same  as  that  obtained  in  section  5.1. 

5.3.   Method  III,  The  Matrix  Solution 

The  matrix  method  which  is  used  with  the  max-flow  problems, 
is  based  on  the  fact  that  any  network  with  n  nodes  can  be 
represented  by  an  equivalent  matrix  of  order  n.   The  basic 
princdple  of  this  method  is  precisely  the  same  as  that  of 
labeling  technique.   The  steps  of  the  method  can  be  summarized 
as  follows: 

Step  1.   Using  the  values  of  excess  capacities  as  the 
elements,  form  an  n  x  n  matrix  where  n  is  the  total  number  of 
nodes  of  the  network  and  allow  two  additional  columns  for  8j  and 
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(1,3) 


(2,3) 


3^ 4   5 

(2,1)  *>>  b>i  (1,3) 


Pig.  21)..   Adjusted  excess  capacities  with 
nodes  labeled. 


3,5=i    g5,3" 


Pig.  25.   Final  network,  no  labeling 
is  possible. 
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Step  2.   Inspect  row  1  for  g-j_ .-  >  0,  write  in  row  j  under 
column  5  these  values  of  g^ ,,    respectively;  and  under  column  r 
write  a  1  which  indicates  that  row  j  is  reached  from  the  source. 
If  the  sink  is  not  one  of  the  j's,  go  to  step  3;  and  if  it  is, 
go  to  step  ;>• 

Step  3.   Consider  the  row  of  lowest  index  from  the  set  of 
rows  labeled  in  step  2.   If  this  is  row  i,  look  for  g. .  >  0  for 
which  row  j  has  not  been  labeled.   Then  set  6  j  =  min(gj_j,  6^) 
and  r-  =  i.   Repeat  the  process  until  all  the  rows  labeled  at 
step  2  are  exhausted.   If  the  n^*1  row  is  labeled,  go  to  step  5> 
if  not,  go  to  step  i].. 

Step  Ij..   Repeat  step  3  for  the  set  of  rows  labeled  in 
step  3-   In  a  finite  number  of  steps  one  of  the  following 
stages  is  reached: 

(1)  So  additional  rows  can  be  labeled,  and  the  n 
row  is  unlabeled. 

(2)  The  n*11  row  is  labeled. 

If  state  (1)  is  reached,  the  existing  flow  is  maximal.   If 
state  (2)  is  reached,  the  existing  flow  can  be  increased;  go 
to  step  S-- 

Step  5-   If  6n  is  the  label  for  the  xfi1   row,  then  deter- 
mine the  path  by  tracing  backwards  on  the  matrix  from  the  value 
of  r  for  the  connected  nodes.   The  increase  in  flow  is  5  . 

Step  6.   Circle  the  elements  on  the  path  in  step  <-,;    note 
diagonally  opposite  elements  and  draw  a  square  around  them. 
Subtract  the  increased  flow,'  5  ,  from  each  circled  element  and 
add  5   to  each  element  with  a  square. 
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Step  ]_.      Repeat  the  process  for  the  new  matrix  formed  at 
step  6. 

To  illustrate  the  above  method,  the  excess  capacities  from 
Pig.  23  are  noted  in  a  matrix  as  shown  in  Table  11.   Then  ac- 
cording to  step  2,  the  positive  elements  in  row  I4.  correspond 
to  columns  2  and  3  with  g-jo  =  3  and  g-,,  =  2,  respectively. 
These  values  are  noted  in  the  6  column  by  rows  2  and  3  and  r  in 
both  cases  is  equal  to  1  since  the  source  is  row  1.   This  ter- 
minates step  2.   Since  5   has  not  been  assigned  in  step  2,  then 
according  to  step  3,  the  row  with  the  lowest  index  from  the 
rows  labeled  in  step  2  is  selected  and  is  the  second  row.   The 
columns  lj.  and  $   can  be  reached  since  the  elements  are  S,o]±   =  1+ 
and  g2cj  =  5,    and  where  the  rows  1+   and  $   have  not  been  labeled. 
Accordingly,  the  rows  l^  and  5  are  labeled  with  6l  =  min(g2j,,  62) 
=  min(I|.,  3)  =3,  r.=  2,  and  5^  =  rain(g2£,  62)  =  mint  5,  3)  =  3, 
r^  =   2,  respectively.   Continuing,  the  next  lowest  row  index 
labeled  in  step  2  is  row  3-      But  there  is  no  column  which  has 
a  row  which  is  unlabeled.   Thus  step  3  is  terminated.   It  is 
noted  that  the  'n"1  row  is  not  labeled;  therefore,  according  10 
step  I4.,  the  first  row  which  has  been  labeled  at  step  3  and  is 
to  be  considered  is  row  I4..   The  column  for  which  the  row  is 
unlabeled  is  column  6.   Now  row  6  is  labeled  by  5^  =  mintgh^,  5r) 
=  mint 6,  3)  =3,  rg,  =  1^.   Finally,  5n=£  =  3;    therefore  flow  can 
be  increased  by  3.   (Refer  to  Table  11.)   Then  according  to 
step  5,  the  path  is  traced  as  follows:   r^  =  i(.,  go  to  row  1± 
and  circle  element  g.  /  =  6.   Then  note  that  r,     =  2  and  go  to 
row  2  and  circle  the  element  g2j.  =  i|..   Note  that  ro  -  1  and  go 
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Table 

11. 

Matrix  of  excess  ca 
( first  iteration) . 
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3 

6 

0 

0 

0 

E 
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k 

Ta 
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1 

- 
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0 
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3 

- 

1 
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5 

0 
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3 

3 
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- 

© 

l 

0 

2    1 

2 

1*. 
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k 
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l 
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3 

5 

0 

k 

0 

0 
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3 
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3 

6 
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0 

0 

12 

0 

- 

2     k 

k 
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1 

and  circle  the 

element 

E12 

=  3. 

Thus  the 

path  is 

l-*2, 

2- 

k,  k 
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and 

is  given  by  the 

circled  e 
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fer  to 

Ts 
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marked  'by 

squares . 

Continuir 

g,  56  is 
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d 

from 

each. 

circled  e] 

ement  and 
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to  the 
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elements  in  squares,  thus  forming  the  new  matrix  which  is 
illustrated  in  Table  12.   The  process  is  repeated  to  determine 
if  the  flow  can  be  increased,  and  continues  until  no  more  rows 
can  be  labeled  and  the  nth  row  remains  unlabeled.   This  then 
indicates  that  the  flow  has  reached  its  maximal  value.   The 
second  iteration  of  the  example  is  illustrated  in  Table  12  and 
the  final  is  illustrated  in  Table  13-   The  sum  of  the  elements 
in  the  last  row  of  the  matrix  gives  the  value  of  maximal  flow. 
Prom  Table  13  it  is  noted  that  the  maximal  flow  for  the  above 
example  is  5  and  is  the  same  as  the  one  obtained  for  this  ex- 
ample by  other  methods. 

Table  13.   Final  matrix  (maximal  flow  =  $) . 


Nodes 

1 

2 

7 

14. 

5 
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Labels 

Step 
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- 
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0 
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1 

1 

5 

0 

3 

2 

1 

- 

0 

1 

0 

k 
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k 

k 
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1 

5 

0 

k 
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0 

- 

3 

6 
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0 

0 

$ 

0 

- 
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5.\+.      Discussion  of  Methods  of  Section  5 

All  the  above  methods  used  to  solve  the  maximal  flow  prob- 
lem are  the  application  of  max-flow  min-cut  theorem,  which 
apply  to  steady-state  flow.   When  the  network  is  large  the 
matrix  solution  eliminates  the  need  of  drawing  the  network  for 
each  iteration  which  saves  a  lot  of  time.   The  example  problem 
considered  consists  of  only  one  source  and  one  sink  but  it  is 
possible  to  have  many  sources  and  sinks.   However,  such 
problem  can  be  reduced  to  the  case  of  one  source  and  one  sink 
(Ford,  Pulkerson  (6),  Hadley  (8)). 
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SUMMARY  AND  CONCLUSIONS 

The  purpose  of  this  report  is  to  review  the  technique 
for  solving  transportation,  transshipment,  and  maximal  flow 
problems.   Some  definitions  and  theorems  concerning  the  theory 
of  network  and  maximal  flow  problems  are  presented.   In  sec- 
tion 3,  several  methods  are  discussed  for  solving  a  simple 
transportation  problem.   These  are: 

1.  Simplex  algorithm 

2.  Transportation  algorithm 

3.  Discrete  maximum  principle 
1]..  Method  of  networks. 

The  methods  discussed  in  this  report  are  only  applicable  to 
problems  with  linear  cost  functions  except  the  discrete  maximum 
principle.   It  is  noted  that  transportation  problems  involve 
shipping  commodities  from  supply  points  to  demand  points  such 
that  the  shipping  costs  are  minimized.   The  transshipment 
problem  contrasts  with  this  in  that  intermediate  points  may  be 
utilized  to  minimize  the  shipping  costs;  consequently  the 
transshipment  problem  is  a  special  case  of  the  transportation 
problem. 

The  method  of  chains  of  positive  arc  capacities,  the 
labeling  technique,  and  the  matrix  method,  which  are  all  re- 
finements of  the  maximal  flow  theory  of  networks,  are  reviewed 
in  section  $,    as  techniques  for  solving  maximal  flow  problems. 
These  problems  are  characterized  by  maximizing  the  flow  from 
a  source  to  a  destination  through  intermediate  arcs  which  may 
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have  limited  capacities. 

In  section  3,  an  example  of  a  transportation  problem  with, 

two  sources  and  three  destinations  was  solved  by  the  various 

methods.   The  number  of  iterations  required  by  each  method  is: 

Method  Iterations 
Simplex  algorithm  l± 

Transportation  algorithm  1 

Discrete  maximum  principle  1 

Method  of  networks  1 

Prom  these  results  and  the  ease  of  programming  the  algorithms 
for  the  computers,  the  transportation  algorithm  appears  to  be 
better. 

In  section  [\.,    an  example  of  a  transshipment  problem  with 
two  sources,  three  destinations,  and  three  intermediate  points 
was  formulated  and  solved  as  a  linear  programming  problem,  and 
was  solved  using  the  transshipment  algorithm  and  the  method  of 
networks.   It  is  noted  from  the  results  that  the  method  of  net- 
works is  fastest  when  the  problem  is  solved  manually.   But  to 
work  out  the  problem  on  computers  the  transshipment  algorithm 
is  the  best. 

In  section  5,  an  example  of  maximal  flow  problem  was  solved 
using  the  method  of  chains  of  positive  arc  capacities,  the 
labeling  technique,  and  the  matrix  method,  which  are  basically 
the  application  of  the  max-flow  rain-out  theorem.   It  seems 
that  the  matrix  method  in  this  case  is  the  fastest  one  as  it 
eliminates  the  need  to  draw  the  network  at  each  iteration. 
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The  purpose  of  this  report  is  to  review  the  techniques 
for  solving  transportation,  transshipment,  and  maximal  flow 
problems.   The  methods  which  are  discussed  for  solving  the 
transportation  and  transshipment  problems  are: 

1.  Simplex  algorithm 

2.  Transportation  algorithm 

3.  Discrete  maximum  principle 
i|..  Method  of  networks. 

It  is  noted  that  a  transportation  problem  involves  ship- 
ping commodities  from  supply  points  to  demand  points  such  that 
the  overall  cost  is  minimized.   This  contrasts  with  the  trans- 
shipment problem,  in.  that  intermediate  points  may  be  utilized 
to  minimize  the  shipping  cost. 

The  method  of  chains  of  positive  arc  capacities,  the 
labeling  technique  and  the  matrix  method,  which  are  all  ver- 
sions of  the  method  of  networks,  are  reviewed  as  techniques  for 
solving  maximal  flow  problems.   The  objective  of  this  type  of 
problems  is  to  maximize  the  flow  from  a  source  to  a  destination 
through  intermediate  arcs  which  may  have  limited  capacities. 
Examples  of  transportation,  transshipment,  and  maximal  flow 
problems  are  solved  using  the  various  techniques. 


